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ABSTRACT 


The sound distribution everywhere within a wedge-shaped fluid overlying either 
a slow or a fast bottom has been studied in this research. Collecting all the results 
of the previous works in this area and overviewing them has been the primary 
purpose. All the cases reported earlier have been studied, and some new ones have 
been added. The variation of the transition point distance with the shore distance 
as a variable was observed. In addition the isopressure patterns were verified by 


calculating the pressure amplitudes in axial direction.. 
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I. INFRODUCTION 


sound propagation in an ocean with a sloping bottom has received considerable 
attention during recent vears. Both tHeoretIclaNs and expdrimentalists have been 
involved with this subject [Ref 1-19]. They have created a number of acoustic 
models to predict the sound field within a wedge shaped fluid overlving a penetrable 
or a rigid fluid bottom. 

The acoustic field in a wedge shaped shallow water duct with ideal boundary 
conditions has been studied by Bradley and Hudimac/Ref. 1]. They have analyzed 
the case of an isospeed duct with one pressure release and one rigid surface. The 
theoretical analysis has been carried out in both image theory and normal mode 
theory. The relationship between the two has been explored for the case of a point 
source. The image sum was transformed via the Poisson transformation formula to 
the normal mode sum. The average acoustic field was then compared for the two 
approaches with good results. Also the two theoretical approaches were compared 
with experimental data from the west coast of Florida with fair results. 

Kuznetsov [Ref 2] used the method of images to investigate the high—frequency 
field in a wedge. He showed that the method of images is better suited to the near 
field where, in addition to the normal modes of the discrete spectrum. the 
continuous spectrum plavs a significant role. At distances from the source greater 
than ten times the laver thickness. the normal mode method works better. ‘These 
results are based on the idealized model of a wave channel as a layer of water having 


plane boundaries over a fluid half space. The sound field in a wedge was described 


as the sum of the image—source fields taking account of the coefficient of all 
reflections. A point source was used in a wedge with perfectly reflecting boundaries. 

Kuznetsov in another study [Ref 3], has described the displacement of the 
normal modes which takes place in the vertical plane and causes each norma! mode 
to emerge from the wedge into the half space at a certain distance. specific to the 
normal mode, from the vertex line of the wedge. 

In 1975 Graves, Nagl, and Uberall [Ref 6] used the adiabatic range variation 
method. derived by Pierce and Milder [Ref 4,5], to perform an approximate 
separation into normal modes for the problem of an under—ocean channel with: 
gradual range dependence of medium and boundaries. Their technique was 
illustrated by application to the isovelocity ocean wedge with rigid ocean floor. and 
compared with the exact solution of Bradley and Hudimac [Ref 1]. A_ good 
agreement was obtained for moderately large wedge angles. even when mode 
coupling was neglected. 

In 1978. Coppens. Sanders. Ioannou and Kawamura [Ref 7]. using computer 
models. predicted the acoustic pressure amplitude and phase in the upslope direction 
at the bottom of a wedge shaped fluid laver overlaying a fast fluid. The slopes of 


the wedges studied were about DES ee DP. 


One of the models had infinite source 
distance while the other had the source at finite distance. The program outputs 
were compared with the experimental results using as variables the wave number 
and density. There were significant differences between the simple model and the 
experimental results because of the experimental! limitations. 


The results of shallow water acoustic experiments performed off the coast of 


Corpus Cristi, Texas. have been compared with theoretical computations by Rubano 


(1980) [Ref 8]. A four—layver. fluid, normal—mode model was used to predict tli 
group velocity dispersion and spatial amplitude distributions of the first and second 
modes. The results show good agreement between the measured and theoretical 
values. 

In the same vear (1980), comparison between experiment and theory was made 
by Jensen and Kuperman [Ref 9]. They used the parabolic equation model to study 
the modal cutoff during upslope propagation in a wedge-shaped ocean. They found 
that the cutoff is not abrupt but takes place over a finite distance which essentially 
provides an aperture for radiation of a beam into the bottom. Also the Gaussian 
representation of a point source in the PE model has given very good results in the 
presence of a bottom. 

The modal cutoff during upslope propagation has also been studied by Pierce 

(1982)/Ref 10]. Since the adiabatic mode theory for upslope propagation in shallow 
water breaks down when the depth decreases to a critical value. he has used 
matched asymptotic expansion techniques to vield the acoustic pressure in the 
transition case. The derived characteristic critical-depth transition function give 
results which are in accord with computations reported by Jensen and Kuperman 
[Ref 9]. 
The same subject, the modal cutoff during upslope propagation. was studied by 
Arnold and Felsen using the image theory as a convenient starting point for 
collective treatment of ray fields and their conversion into local modes [Ref 11]. 
The numerical comparisons with the results of Jensen and Kuperman [Ref 9] showed 
good agreement. 

The sound field in an absorbing fluid substrate underlying a wedge shaped fluid 


with higher sound speed has been studied by Coppens, Sanders, and Humphries in 


1984 [Ref 12]. A Green's function extension of the image theory give the 
opportunity of doing all the calculations on a computer. Two different 
approximation methods were used. The end—point approximation gave solutions 
valid at great distances. and for absorptions representative of sedimentary 
materials. The saddle—point approximation has given good results under the apex. 
with no limit in the absorption. 

In 1984 Baek [Ref 13] predicted the pressure amplitude and phase, everywhere 
within the wedge, in upslope direction, having a fast bottom. Three cases were 
studied. pressure release bottom, rigid bottom, and penetrable bottom. 

In the same vear (1984) LeSesne [Ref 14] studied and compared two computer 
programs using the method of images to determine the pressure and phase 
distribution. It was concluded that an infinite source is not approached until the 
source exceeds six hundred dump distances from the apex. Al! measurements were 
taken directly upslope with source distance equal to forty dump distances and the 
receiver on axis with the source. Comparison between the two programs show 
similar results. 

One year later, in 1985 Borchardt [Ref 15] experimentally measured the sound 
pressure field everywhere within isospeed water overlying a sloping absorbing 
bottom. The experimental data were compared with the program developed by 
LeSesne [Ref 14]. The results showed that the program successfully predicts the 
on—axis upslope sound profile in a wedge-shaped medium. Particularly, 
measurements showed good agreement for R2< 10, 6 = 9.5°, P1/ Po = |) ae 
Y0<50. All the symbol definitions are given in the List of Symbols table. 

A shallow water wedge with a slope up to 9° has been experimentally studied 


by Tindle. Hobaek, and Muir in 1986 [Ref 16]. Thev presented experimental results 


for downslope propagation and they defined the wedge modes as propagating normal 
modes with wavefronts which are curved into arcs of circles centered on the wedge 
apex. 

In 1987 Kaswandi [Ref 17], studied the slow bottom case. The outputs were 
taken for the on—axis case with wedge angles of 6°, 10°, and 15°. Three kinds of 
curves were observed. depending on the source distance. 

At about the same time, a program with 3—D capabilities was used by Li [Ref 
18] to predict the acoustic pressure distribution along the interface between a 
wedge-shaped fluid laver overlying fast or slow bottom. It was concluded that the 
major characteristics of the output graphs were insensitive to minor variations of 
the bottom. 

In 1988 Doolittle and Tolstoy presented experimental data obtained in East 
Australian Continental Slope consistent with the theory of energetic horizontal 
refraction due to multiple reflections from a sloping bottom [Ref 19]. When viewed 
from above, the path of the ray appears to curve. As a result,the exact position of 
the source can be obtained from the details of this curvature. 

The purpose of this research is to take all the outputs provided by the theses of 
Kaswandi and Li and trv to build a physical understanding of the acoustic field from 
them. Also an attempt was made to obtain detailed comparisons between the 
DSLOW and 3LS4 programs for the same inputs,and to further investigate the 


transition point. 


I]. THEORY 


The theoretical background of this research has been presented by Coppens. 
Sanders, Joannnou and Kawamura [Ref 8]. What follows, in this chapter, is a 
summary of the most important points of the image theory and its apphcation to 


the wedge-shaped ocean with penetrable bottom. 


A. METHOD OF IMAGES 

The assumptions that (1) the speed of sound is the same everywhere within the 
wedge and in the penetrable bottom and (2) the upper surface is a pressure release 
surface makes the method of images an appropriate approach for understanding the 
sound field in a wedge-shaped fluid. The sound paths are straight lines and the 
plane wave Ravleigh reflection coefficients. discussed by Kinsler. Frey. Coppens 
and Sanders [Ref 14], is a good approximation if the sound source is not too close to 
the bottom. 

The geometry of the model used to predict the pressure amplitude and phase 
evervwhere within the wedge is shown in Figure 1. For the wedge-shaped duct. 
cvlindrical coordinates are used with the shore line as the axis (Figure 2). The 
images lie on a circle centered on the shore line. 

In Figures 3 and 4. the relationship between the receiver and the nth image is 


shown. The angle 6, of the nth image, measured from the bottom, is given by: 


$, = B-(n-1)+G, for n odd 


é,= B-n—-G for n even 


According to Sommerfeld's concept of an extended Riemann surface. thie 
extension of the @ coordinate makes the number of images infinite. The tota! 
pressure and phase at any field point within the wedge is found by coherent 
summation of all spherical waves radiated from the infinite number of images. As 
the number of reflections associated with a given image increases, the pressure 


contribution of this image decreases. The infinite summation is approximated by 


taking the sum over N images where 
ele Gein liSOs) 2) 
Wedge angles of the form z/n have been used, where n is an integer. because 


diffraction terms were eliminated [Ref 2]. 


The distance Ry between the receiver and the nth image is given by: 


Ry-+ Ro? — 2-RyRo-cos(@—D) + Yo" 





for the upper group of images, and 


Rn -cos(@n+D) + Yo" 





for the lower group images. 


The angle of incidence Opn of the nth image on the m"® plane ( described by the 


angle m: ) is 


=—| 


([Ry-sin(@:— m3)+ Roe-sin(m3 — D)] 
SHO — 
Ry 


The reflection coefficients R(@,) for the upper family of images and R(@,') for 


the lower family of images for a plane wave are 
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" 1 = (c1/c2)*cos?( Onn) 
nm = 
sin( @nm) 
and 


y. 1 — (c;/c2)?cos?(On'n) 


nm = 
sin(On'm) 


The resulting complex pressure from the upper family of images is 


and the resulting complex pressure from the lower family of images is 


ns 


eS = exp(—jkRq')- (1) NM+)/2] . Tr R(a,',) 
n=] Rn on) 


The sound field radiated by each image has been multiplied by the plane wave 
reflection coefficients corresponding to reflections encountered by the wave as it goes 
from the source to the field point. 


The total complex pressure distribution P is the sum of the Pu and P): 


P= Pu + Pj 


B. NORMALIZATION 

All distances are normalized in terms of the characteristic distance X.. This is 
the distance measured from the apex along the wedge interface to the point at which 
the lowest mode attains cut off as shown in figure 5 [Ref 20]. 

The characteristic distance is defined as the distance from the shore at which 


the depth H is given by 





to 
— 
| 
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For a fast bottom. the critical angle @ is 


CLs 
<n cos bc 
since 
i 
Soe — 

2sin Oe 

tang = H 

x 

t 

kX = 


2 sin|@,}| tang 


The above formula gives the derivation of the characteristic distance X in the 
case of a fast bottom (C) < cy). When slow bottom case is studied (Cj > Cp aa 
characteristic distance can not be defined as above. A useful characteristic distance 


for slow bottom can be defined by inverting ¢1/c2 


kX = i 
2 tanlarcos(C,/C, )] tanZ 


The advantage of working with normalized distances is that the analvsis is 
frequency independent. 

A convenient normalization for the pressure amplitude is accomplished by 
multiplying the pressure amplitude at one meter from the source by the source 
distance (Rl). This simplifies the interpretation considerably, especially in the case 


of large source—apex distance compared to X. 
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In the course of this research. a computer program has been created. the DEML 
program. This computer program provides an additional normalization. |i 
normalizes all the pressure amplitudes to the first pressure measurement taken down 
from the surface. This slope normalization makes the shapes of the curves 


comparable, even if thev have been obtained at different distances. 
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Il. COMPUTER PROGRAMS 


Two programs have been used for the purpose of this research. They are 
improved and more documented versions of the programs discussed by haswandi 
[Ref 17] and Li [Ref 18]. They use the same model but have different output 


displays. 


A. PROGRAM DEME 

In the beginning of the study reported in this thesis, the DSLOW program was 
obtained [Ref 17]. DSLOW is itself an improved version of the WEDGE and 
XSLOPE programs developed by Kawamura and Ioannou [Ref 7]. and LeSesne [Ref 
14] respectively. 

DSLOW was developed to run on the IBM 3033 main frame of the Naval 
Postgraduate School. It uses the method of images to predict the pressure amplitude 
and phase within the wedge fluid overlying a slow bottom. 

DEME. created from DSLOW, uses the same model to predict the pressure 
amplitude, and phase within the wedge fluid overlving a slow or a fast bottom. 

The output is given in a table format or in a graph as shown in the Appendices. 
The DISSPLA graphics package and the EASYPLOT program, both of them 
available on the IBM 3033. have been used for the graphical outputs. 

During this research, other versions of the original DEME program have been 
created with different outputs. All of them have been stored under Professor 


Coppens’ account in the IBM 3033 in the Naval Postgraduate School. 


The most important of these programs are: 

1. Program 5 provides the pressure amplitude and phase at one specific poiut ir 
the field. It must be used with the "Record on" command of the IBM 3033 for 
sending the results in a separate file. 

2. Program SDF is a version of the DEME program which takes the input data 
from an input file and transfers the output to another file. With this technique the 
EASYPLOT program can be used to produce multiple graphs. 

3. Program PREC uses double precision to reduce round—off errors. The outputs 
of PREC can be plotted with any other graphics software except DISSPLA. The 
DISSPLA package has some difficulty when applied to the output of double 
precision fortran programs. 


4. Program SLAI generates graphs without the slope normalization. 


B. PROGRAM 3LS4 

This program gives alternate forms of data presentation. It uses the IBM 3033 
and the DISSPLA package to generate a 3—D representation of the pressure 
amplitude distribution on the bottom or along any plane (defined by the D) within 
the wedge. It also provides a contour presentation of the same data. 

Because of the large volume of numerical output. 3LS4 is ineffective for 
recognizing details. 

An important limitation of the program is that the array containing the data 
must be initially dimensionalised to the exact size of the data file array. This 


information is normally passed from the data file during the program. 
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IV. DISCUSSION OF THE RESULTS 


The DEME and 3L54 programs have been created using the same model. 
Although they provide different output graphs they must agree at every point. The 
3LS4 provides the pressure amplitude and phase on a plane defined by the receiver 
angle (D). The DEME program provides the pressure distribution and phase along 
a vertical line anywhere within the wedge. The programs were checked for both 
slow and fast bottoms. For the same input data, thev gave the same outputs. as 
shown in Tables 1 and 2. 

For a slow bottom (c; > Ca), three different types of curves of pressure amplitude as 
a function of depth have been observed. The same types of the curves have been 
noticed by Kaswandi [Ref 17] and are shown in Figure 6. 

When the source distance (R1) is much smaller than the characteristic distance 
(X). the sound pressure is linearly dependent on receiver angle with the surface 
pressure equal to zero and the bottom pressure a maximum (Rl = 0.1) When RI 
approximates X , tlle curve develops considerable curvature (R1 = 0.7). 

When the source distance (R1) is slightly different from the characteristic 
distance X (R1 = 1.1), then a maximum and a minimum appear between the 
surface and the bottom. 

When the source distance (R1) is greater than the characteristic distance (X) 
(Rl = 4). then the sound pressure has a maximum between the surface and the 
bottom. No minimum appears in this kind of the curves. When the value of R1 is 
much greater than X , several minima and maxima appear in the output curve 


(= 2). 
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Figure 7 shows the development of the second type of curve (R1 = 1.0) as the 
receiver distance is varied. 

E@ueume saine input data (B= 10, G = 5, p/P, - 0.9), the fast bottom gave 
entirely different curves as shown in Figure 8. When the source distance (R1) is 
less than three times the characteristic distance (X), only the first mode is received 
everywhere within the wedge (R1 = 0.3, 1.1, 1.3). When R1 > 3. higher modes are 


excited and the curves are the superposition of several modes (R1 = 4). 


A. TRANSITION POINT 

For the slow bottom, when the source distance (R1) is about the same as the 
characteristic distance (X), the curve has a maximum and a minimum. The 
minimum is close to the bottom. A "transition point" has been defined bv 
Kaswandi [Ref 17] as the receiver distance, for which the pressure above a minimum 
extrapolates to zero on the bottom. 

Three examples are shown in the Figure 9 (R2 = 4.6, 6.4, 32). In all of them 
the pressure above the minimum can be extrapolated to zero pressure point at the 
bottom. The transition point phenomenon happens for source distances 1<R1<1.5. 

For 1 < Rl < 1.5 the shape of the pressure curves as the receiver distance 
(R2) is varied can be explained as follows: as the receiver is moved from (R2 = 1.5) 
to (R2 = 30) the minimum pressure amplitude decreases in magnitude and the 
first transition point is observed (Fig 7). With further increasing of R2. the 
minimum pressure amplitude of the curves decreases to a minimum value, and then 
increases arriving at a second transition point. For the receiver distances between 
the first and the second transition point, the extrapolation to zero pressure was 


observed to be above the bottom (Figure 7). With further increasing of R2. a third 


transition point tries to form but. since the R2 is so large round—off error appears 
in the program output curve and the curve characteristics can not be easily 
recognized. In the figure 10, the locus of the two detectable transition points are 
shown as a function of R2 and the shore distance (YO). Although it is not indicated 
in the figure 10 the points lie on circular arcs centered at the source. This may 
indicates that the transition point obevs the cylindrical spreading law. 

The locus of the first and second transition points form curves which can be 
approximated by arcs of circles centered at the source (figure 10) 

The transition point as a function of the source angle (G) was studied using the 
output data shown in Table 3. As the source angle changes, the transition point 
moves to a different receiver distance. Two different wedge angles were observed. 
The greater the source angle (G) . the smaller the receiver distance where the 
transition point occurs. 

A possible explanation is that, by increasing the G the source is placed far from 
the angle value (D), where the minimum happens, so the absolute pressure values 
become smaller as the source angle increases. In figures 10a and 10b the pressure 
distribution of the bottom is shown for G = 2.5°. In figures 10c and 10d the same 


case is represented with G = 7.5. 


These figures are generated with the program 
3LS4 [Ref 18]. 

In the fast bottom case no transition point was observed. In Figure 11 the fast 
bottom case is shown when the shore distance (YO) varies. This is the case of 
downslope propagation in a wedge-shaped fluid overlying a fast bottom. At the 


receiver's position (R2 = 32), as the shore distance is increasing ( 0 < YO < 80). 


only a curve like the first mode appears and no minimum is observed. 
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B. MORE CASES 

For the fast bottom, when the source distance (R1) is bigger than 3. the higher 
modes start to appear. In Figure 12, Rl =6. As explained by Kaswandi [Ref 17| 
when R1 = 6 the first three modes are present. By increasing the source angle (G) 
the output curve changes because at some source positions some modes are reduced 
and some are amplified. 

some problems appear when the source is very close to the apex (R1 =0.1). as 
in Figure 13. The program seems to fail close to the bottom. It may be a round—off 
error of the computer because the absolute pressure values are very small. The 
round—off error is clearly shown in Figure 14, where it starts to appear when the 
source—receiver distance becomes too far. An attempt to eliminate this problem 
was made in the PREC program which was a double precision version of the DEME 
program. Also another try was made to run the DEME program in a 386 IBM 
compatible personal computer with 16—bit memory. Both of the tries gave 


smoother curves but they did not solve the problem. 


C. SHORE DISTANCE AS A VARIABLE 

Until now most studies were for on axis receiver positions (YO = 0). In the 
following discussion the shore distance (YO) is the variable. In slow bottom three 
different source distances (R1) were studied (R1 = 0.5, Rl = 1 and Rl = 5). 

The first case is shown in Figure 15, where the source is close to the apex (R1= 
0.5). At receiver position R2 = 4. only the linear pressure curve is observed at all 
values of YO (0 < YO < 60).The effect of larger values of source and receiver 
distance (Rl= 1, R2= 6) is shown in Figure 16. As YO increases, the curve forms a 


transition point at YO = 3 and then becomes a straight line at YO = 10. 


in 


Figure 17 illustrates the case when the source distance is greater than the 
characteristic distance (Rl = 5). Curves for four different values of YO are shown. 

The same cases were studied also for the fast bottom. The outputs are shown in 
Figure 18 when the source distance is close to the apex(R1 = 0.5). in Figure 19 
where Rl = 1, and in Figure 20 where the Rl = 5. What is clear to be observed in 
these cases is the uncertainties. which are present in great distances. Also when 


R1> 3 the higher modes start to appear. 


D. ANGULAR OUTPUTS 

Examination of the contour plots revealed a general hyperbolic or elliptical 
pattern, therefore another trv was made to understand the physical structure of the 
acoustic field by searching at some particular angle off the axis ( arctan[Y0/R2] ). 
As shown in Figure 21 output data were calculated every 10° off the axis. The 
acoustic field was studied for two different source positions, for a small source 
distance(R1=0.5) and for a source distance greater than the characteristic distance 
EOL SIMOHe 

The bottom pressure distribution is shown in Figures 22 through 25. The slow 
bottom as well as the fast bottom has been studied. From the contour plots 
(Figures 23 and 25) it would appear that large values of the axis angle were more 
interesting. because the isopressure curves seem to be asymptotic to the axis angles. 

In the Figures 26 through 28 , the pressure distribution along three different 
angles are shown for the slow bottom. The curves are unnormalized. 

The fast bottom graphs for the same inputs are shown in Figures 29 through 34. 
Many different axis were studied from 0° to 90°. All the curves from 0° up to 70° 


seem to be the superposition of the first five modes (Figure 29, 30,and 31). From 
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75” axis angle and for big receiver distances (R2) the graphs begin to have only thie 
first mode. That means that the higher modes have been cut off. At 80° all the 
Outputs appear to be in the area where only the first mode is received. The R2 in 
all the graphs was varied between 10 and 60. 

The slow bottom curves for any angle off the axis were observed to have 
uniform vertical pressure profiles as shown in Figures 26 through 28. The fast 
bottom curves for angles up to §0° were observed to have complicated profiles as 


shown in Figures 29 through 32. 


D. THE OTHER PARAMETERS 
Variations of the other parameters such as the velocity ratio c;/co and the 
density ratio p, /fP. were studied by Kaswandi [Ref 17] for on axis measurements. 


The results showed no strong variation of parameters studied. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


In this research. the 3—D propagation of sound has been studied using the image 
theorv. Two computer programs were used to provide the pressure amplitude 
distribution and phase everywhere within a wedge-shaped fluid overlving a fast or a 
slow bottom. 

The DEME and the 3LS4 programs were checked for the same input data and 
thev provided the same outputs for the same acoustic field points. After this check. 
it was concluded that the programs were consistent with each other. 

The transition point, present only over the slow bottom, was studied at various 
places inside the wedge. It was found that all the curves could be approximated by 
concentric circular arcs. having the source at the center. 

As the source angle (G) increases, it was concluded that the transition point 
distance decreases as it is shown in Table 3. 

When the source position is 1.1< Rl < 1.5 then two transition point curves 
were formed. Sometimes it was possible to have also a third one, if the computer 
could work with very small numbers (16 or 32 bit Bus). 

The off axis data were taken along many angular axes. The slow and the fast 
bottom were studied for the same angles of axis. The slow bottom data showed very 
small variation of the pressure amplitude values as the receiver distance was 
increased. The fast bottom curves showed very clearly the wave guide cutoff 
phenomenon. In all the axis angles up to 80° the curves were observed to be 
irregular with many variations and unique for each R2. At the angle of 80° or 


more, nO major variations were observed as the R2 was increased. This observation 


was found to be consistent with the contour plot output data taken from the 3LS! 
program. 

This research was another phenomenological] theoretical step toward elucidating 
features of the wedge problem. Some experimental! data must be compared with tlie 
two programs DEME and 3L54 to have some confirmation of the programs 
capabilities. The most interesting case seemed to be the slow bottom Case. 

Current thoughts [Ref 23] suggest that these features can be studied from 
another point of view based on a collection of dipoles. The source and its first 
image form a dipole as do other associated pairs of images (Figure 1). If al] these 
pairs of images are considered, then the field can be considered to be formed from 
the near cancellation of equivalent dipole pairs lving at slightly different distances 


from the receiver. 
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APPENDIX A 


TABLE 1 





CONSISTENCY OF THE DEME AND 3LS4 PROGRAMS (SLOW BOTTOM) 
B= 10°,G=5°, YO = 50, Rl = 40, R2 = 10, pl/p2 = 0.9, cl/c2 = 2 


REECE ER OUTCOMING PRESSU Rite AM PIA sis 
ANGLED JON VL 3LS4 

00 0.5073 0.50728 
01 0.8683 0.86825 
02 1.4251 1.42509 
03 Wokelehe 1. 396ou 
Ov 1.8491 1.84912 


nt) 0.0000 0.00000 


ie) 
Le) 


TABLE 2 


CONSISTENGY_OF THE DEME AND 3L54 PROGRAMS (EAST BOTTOM 


B= 10°.G =5°, R1 = 40. R2 = 35, p;/P) = 0.9, C,/C. = 0.9 


MeeeANGLE SHORE DIST. 


D Y0 
00 00 
00 30 
00 00 


PRESSURESANMPLI IDE OlUEPUT 


DEME 3LS4 
9.5446 9.54476 
2.0 tO 2.57740 
3.2210 3.22089 


TABLE 3 
TRANSITION POINT AS 45RUNCTION OF SOMRGE A NGG) 


P, [Po = 0.8, C,/Cy = 1.1 


WEDGE SOURCE SOURCE RECEIVER 
ANGLE (B) ANGLE (G) DIST.(R1) DIST.(R2) 
10° Da 1.0 33.0 

ie) 5.0° 1.0 24.0 

10° 7.50 1.0 22.0 

06° 1.5° 1.3 5.9 

06° fo 1.3 5.7 

06° 4.50 1.3 5.4 
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Figure 10A. Pressure Distribution along the Bottom (G = 2.5) 
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Figure 10C. Pressure Distribution along the Bottom (G = 7.5) 
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APPENDIX C 
DEME PROGRAM 


—_—_ oO eee 
—_—_— om iT i meee 


This program calculates the pressure amplitude and phase everywhere withis 
a wedge shape fluid overlying a penetrable bottom. It is a modified version of the 
DSLOW program, which was written by CASWANDI. It has been modified by 
LT.D.PALIATSOS on March 89. This program is working for fast or slow bottoms. 
The output data are given by a table or by a graph using Disses: 


* KK x OK *K x KK x KK 


PSE peel tt MgNiS1S2,N1,J,K 


REAL*4 BCC,C2 D,D1,.D2,G, Pl PuPvOlmninn? 1. 
* T4,T6,W0,W1,Y0, Y1,Y2,Z1,Z2,Z3,24,25.26. 
; T1(80),R8(80).R9(80),S(30).C(30),E(30), 


F(30),Y.Z.R3,AL,PZ(30),DZ(30),V,DX.XP,DD,PN(30) 
REAL‘ TQQ,TQQ1,TQQ2,TQQ3 
= ACOS(—1.0D00 


Af. Reciareh WT Sage 2K OK OK OK 2 3k 2K OK OK 2 OK ok ok OK ok KK OK ok KK OK OK KK OK OK KK ok KK De a ok OK OK OK OKO Oe 


feet Neu PARAMETERS 

ti ees EEE RAE EE A RAE Ee Oe See ee Ee = 
B = WEDGE ANGLE (DEG) 

G—SaeURer ANGLE (DEG) 

D = RECEIVER ANGLE (DEG) 

Nl= # OF IMAGE POINTS 

R= OURere DISTANCE (IN DUMP DISTANCES) 
R2— RECEIVER DISTANCE (IN DUMP DISTANCES) 
YO = APEX DISTANCE (IN DUMP DISTANCES) 

D1 = RHO 1/RHO2 

CC] Cae-2 

AL = ALPHA/K2 

heer aechiV ER POSITIONS 


a ee GS ee EG KR RR a eS aS ee 


Cre! Ut sre CT SRECA TIONS 


a ee ERE ET TREK R EE Ke EX ATE Erg ME P= 


CC) Crean) CO) a Ga 


DO 90 P=1 4 
MOU NUS GIVE THE INFORMA TION ABOU AMESBOTTOM, THAT 


Mino y OU MUST CHOOSE SLOWIOR FAST BOTTOM BY GI¥ING 
Seen, ALUE TO THE SPEED RATIO CC =Cl1/C2. 


CC=GiWe@? 


PIO Oe) 2 fe 


(1 He EO 2 ER ee eee eee 


MAIN PROGRAM 


(UA a 2 a REO OK cE ES A Te ene 


Nl = Pa eae 


T6 = 180./PI 

B = B/T6 

G=G/T6 

CG) = Ceo 

D2 = (YO*Y0)+(R1*R1)+(R2*R2) 
R3 = 2.*R1*R2 


C 4 = PI/(2*TAN(ACOS(1/CC))*TAN(B)) 
TQQ = TAN(B) 


c 
C DECISION ABOUT SLOW OR FAST BOTTOM 
‘© 
IF (CC.LT.1) THEN 
TQQ1 =ACOS(CC) 
TOQOQ2 = SIN(TQQ]) 
C 4 = PI/(2*SIN(ACOS(CC))*TAN(B)) ***** FAST CASE 
ELSE 
TQQ1 = ACOS(1/CC) 
TOO? = TANT OO] 
T4 = PI/2*TAN(ACOS(1/CC))*TAN(B) **** SLOW CASE 
ENDIF 
Cc 


TQQ3 = 2.*TQQ2*TQQ 
T4 = gee 

C HDGEVe 1oosie a2) 
D2 = YO"Y0+R1*R1+R2*R2 
Q1 = 1/DSQRT(2.0D00) 

C300 FORMAT(' NOP '.7X,'THETA(N) ',7X."IMGE SR R8 ',7X. 

C * ‘IMGESR RQ’) 

C WRITE (6,300) 

800 FORMAT(5X.'REC.POS'.5X,' REC.ANGLE ',5X,'PRES.AMPLITUDE ',5X. 
* 'PHASE ANGLE '(5X,'NORM.PRESS') 
WRITE (6,800) 

DO 10 M=0,A 
D = M*B/A 
DD=B/A 
V =2*B/10 
Le ey GOTO 110 
IF(D.GE.V) GOTO 120 


oh Seale 


ae 8) 
ares 
et 
=" 
= 
ll 
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C310 FORMAT(3X,1 
WRITE (6,310) N 


C 


=e 

0) 

2) 

HS) 


Ip 
lac dh 
2 


20 CONTINUE 
r= 0.0 

P2 = 0.0 

DO 30 Nee 1,N1 


C 


C 


C400 FORMAT(' NOP ',5X,'I ',5X,'ARSIN(S(I)*T6)',5X,'E(1)' 


C 
C 


Yee 


W1 


T1(N),R8(N).R9(N) 


aa 
Ti ')=(N-1)*B+G 
T1(N)=N*B-G 


D2-R3*COS(T1(N)—D)) 

D2—R3*COS(T1(N)+D)) 
\5X,F5.2,4X,F6.4,6.X,F6.4) 

N 


)™*(INT(N/2)) 
Seo" AL 


11 = INT((N—1)/2) 
pO 40 I=1,1 


S(I) = 
C(I) = SQRT(1 
1)/D1 


ABS(R1* SIN(T1(N)=2*I*B) 


+R2*SIN(2*1*B— aa Re) 


WO = (-C2+(C(I)*C(I))) 
Y = SQRT((WO0*W0)-+(\W1*W1)) 


Z 


=ABS(W0) 
IF(Y.LE.Z) Y = 


—(S(1)*S(1))) 


Z 


Y1 = Q1*SQRT(Y+W0) 
Y2 = —Q1*SQRT(Y-W0) 


Zi) 


= T-Y2 


72 =~-Y1 
Z12= CMPLX(Z1R.Z21) 
73 = ate 


L4 =—Z2 
Z1= T+Y2 
i yi 


Cie +29. 72) 


Z13= CMPLX(Z3R.Z41) 
Z5 = Z1*Z3-22*7A 
Z6 = Z1*244+22*23 
Z14= CMPLX(Z5R.Z61) 


E(I 
Fd 


o ox, 


= 28 
a6 


'F(T)',5X, 


C WRITE (6,400) 


C220 Sae@r NA Tel? 6X12 5X.F6.4,10X,F6. 


C 


* —— ' DSQRT(E(1)*E(I)+F(I)*F(1))',5X,'T6*ATAN(F(I)/E(I))') 
4 


WRITE (6,420) NI, ASIN(S(1) °T 6,E(1), 
* F(] ( 


C —_ SQRT(E(I)*E(1) +F (1)*F 
ANC Omir WE 


| a F6.4,3X,F6.4,10X.F7.4) 
F(I 
)) ATAN(F(I)/E(1)) 


75 = 1 
Z6 = 0 
IF(N.LE.2.00) GOTO 50 
DO 45 ]=1.11 
Z1 = E( 
Z2 = F(I 
73 = 25 
Z4 = 26 


Z5 = Z1*Z3—Z2*Z4 
Live 7)" 74-2 oe 


45 CONTINUE 


50 


Z1 = 25 
Z2 = Z6 

= T4*R8(N) 
Z3 = COS(T) 
Z4 = —SIN(T) 


Z5 = Z1*Z3-22*Z4 
Z6 = 21*Z4+2Z2"2Z3 
Pl = P1+S$2*Z5/R8(N 
P2 = P2+S2*Z6/RS(N 


C500 FORMAT( NO OF I.P 3X, RE(REFL)= ape ee) 


WRITE (6.500) 


C510 FORMAT(3X.12,6X.F6.4.6X.F6.4) 


’ WRITE (6,510) N,S2*Z1.S2*2Z2 
C600 FORMAT(' I'.5X,'ARCSIN(S(I))*T6='/,3X,'E(1)',5X,'F(1)',5X,'EF’, 


5x! ANGE 


WRITE (6,600) 


l1=I1+1 
DO 601=1.11 
S(I) = ARG RAST = -1)°B) 
: G2 SIN ese) Oo 
C(I) = SQRT(1-S(I)*S ()) 


T = S(1)/D1 

WO = —C2+C(I)*C(I) 

Y = SQRT((W0*W0)+(W1*W1)) 
Z = ABS(W0) 

IF(Y.LE.Z) Y = Z 

Y1 = QI*SQRT(Y+W0) 
Y2 = —Q1*SQRT(Y-W0) 
7 = te 

72 =-Y1 

73 = Zt) (ae 77) 
i a a 
Z1 = T+Y2 

Z2 = Y1 

Z5 = Z1™Z3-Z2"Z4 

Z6 = Z1*Z4+Z2"23 


: = 25 

F(I) = Z6 

AN1= ASIN(S(I))*T6 

EF = SQRT(E(1)*E(1)+F(I)*F(I)) 
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AN2= ATAN(F(I)/E(I))*T6 
C610 FORMAT(2X. > 5X. F6.4,5X.F6.4 45x. F6.4,5X,F6.4,5X.F7.4) 
C WRITE (6,610) I,AN1.E(I).F(I),EF 
60 CONTINUE 


Z1=0 
72 = 
73 =0 
74=0 
Z5=1 
76 =0 
DO 801 = 1,1 
11 
2 = FU) 
23 = Z5 
Z4 = 76 


Z3 = Z1*Z3—Z2*Z4 
Z6 = Z1*Z442Z2*Z3 
80 CONTINUE 


Z1 = 25 

Z2 = 26 

T = T4*R9(N) 
Z3 = COS(T) 
Z4 = —-SIN(T) 


ih 1 RE ag Fe 
Z6 = Z1*Z44+2Z2*73 
Pl = P1+$2*Z5/R9(N 
P2 = P2+S2*Z6/R9(N 
C700 FORMAT(' LOWER PATH NO='/3X.'RE(REFL)= '.3X.'IM(REFL)= ') 
C WRITE ( 6,700) 
C710 FORMAT(6X.12,12X.F6.4,12X ,F6.4) 
C WRITE (6.710) N,S2*Z1,S2*Z2 
30 CONTINUE 
k=K4+41 
A aN: 
PZ(K)=SQRT(P1*P1+P2*P2)*R1 
C WRITE (6,810) K.DZ(K),PZ(K).ATAN(P2/P1),PN(K) 
C810 FORMAT(6X,13,11X,F5.2,12X,F7.4,12X,F7.4,10X,F7.4) 
15 CONTINUE 
10 CONTINUE 
DO 31 L=1,29 
XP=PZ(28) 
PN(L)=PZ(L)/XP 
WRITE (6,811) L.DZ(L),P 
811 FORMAT(6X,13,11X.F5. 
31 CONTINUE 
C250 FORMAT(' WEDGE ANGLE = ',F5.2,1X,! SOURCE ANGLE ' 
C * ‘SOURCE DISTANCE='F4.2,1X,'RECEIVER DISTANCE= ', 
C * ‘SHORE DISTANCES ',F4.2) 


Z(L). ATAN(P2/P1),PN(L) 
2,12X F7.4,12X,F7.4,10X,F7.4) 


C WRITE (6.250) B*T6.G*T6.R1,R2.Y0 
Cz0 FORMAT(’  RHO1/RHO2=',F5.2,5X,' C1 /C2='.F5.2.5N. 
AIPHA/K2=',F8.4) 
C WRITE (6.270) D1,CC,AL 
C YO =Y0+2 
C R2= R242 
C R1l=R142 
90 CONTINUE 
B = B*T6 
G = G*T6 
C Sir 
C END 
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C A PROGRAM FOR PLOTTING BY TEK618 OR SHERPA 

(CR ER AI EAA A Ae FEE A RIE RO Ee ena ae 
CALL MEDBUF 
CALE TEKGIS 

C CALL SHER AGSre@n WENN" A’ 3) 

CALL NOBRDR 

CAIIBPAGE(5 5) 

CALL PAGE(S.5,11.) 

CALL HWROT('AUTO’) 

CAE AREA2D(009.,5.) 

CALL AREA2D(09.,7.) 

CALL AREA2D(8.5,6) 

CALL AREA2D(6.0,8.5) 

CALL HEIGHT(.2) 

CALL XNAME('NORM. PRESSURE AMPLITUDES' 24) 

CALL YNAME('RECEIVER ANGLE(DEG)S',19) 

GANLTL YTICKS(S 5 

CALL XTICES 

CALL GRAF(0..2.0 ,10.0.0.,1.0.6.0) 

CALE SRARO. le 0 ,05.0.0.,1.0.06.0) 

CALL Doe 

CALL GRID(2.2) 

CARE HEADIN(REC.ANGLE VS. NORM.PRESSURES',—100,1.8.1) 


Qe Co & 


© 





CALL MESSAG('WEDGE ANGLE =$",100,8.,7.) 
CALL REALNO(B,2,10.2,7.) 

CALL MESSAG('RHO1/RHO2 =$',100,8.,6.5) 
CALL REALNO(D1,2,10.2,6.5) 

CALL MESSAG('C1/C2 =$',100,8.,6.0) 

CALL REALNO(CC,2,10.2,6.) 





CALL MESSAG('SOURCE ANGLE= $',100,8.,5.5) 
CALL REALNO(G,2,10.2,5.5) 

CALL MESSAG('SOURCE DIST.= $',100,8.,5.) 

CALL REALNO(R1,2,10.2,5.) 

CALL MESSAG('REC.DIST. = $',100,8.,4.5) 

CALL REALNO(R2,?2,10.2,4.5) 

CALL MESSAG('SHORE DIST.= $',100,8.,4.) 

CALL REALNO(Y0,2,10.2,4.) 

CALL MESSAG('PRESS.AMPL.(X=1) = $',100,7.0.3.5) 
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CALL REALNO(XP,105,10.2,3.5) 
CALL RESET(‘ALL') 

CALL PARA3 

CALL NOCHEK 

CALL CURVE(PN.DZ,29,1) 
CALL ENDPL(0) 

CALL DONEPL 

STOP 

END 
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